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Abstract 

The linear and renormalized nonlinear analysis of the temporal evolution of drift-type modes 
in plasma flows with strong time-varying velocity shear is developed. Analysis is performed in 
the time domain without spectral decomposition in time and admits time variation of the flow 
velocity with time scales which may be comparable with turbulent time scales. 
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I. INTRODUCTION. 



A great variety of experimental results are now available[l], which confirm a suppression of turbulence 
by E x B velocity shear as a key feature of the regimes of improved confinement [2] of plasma and 
formation of transport barriers in tokamaks. The first remarkable success in theoretical investigation 
of this fascinating and rich topic of fusion research was achieved in Refs.[3, 4], where nonlinear 
enhanced decorrelation of fluctuations by shear flow had been proposed as a universal and robust 
mechanism for turbulence quenching and transport reduction. It was compared quite favorably 
with experimental data for numerous tokamaks and was confirmed by gyrofluid codes[5, 6, 7]. This 
theory was developed with the help of a simple "modal" approach, under an implicit assumption 
that perturbations of fields, densities, currents ets. are "normal eigenmodes" with certain frequencies 
and wave numbers. However shear flow itself affects the structure of the waves and their temporal 
evolution. The shear flow distorts wave pattern, leads to its stretching and to temporal changing 
with time the wave number of any separate wave[8, 9]. Under condition of strong flow shear, which is 
observed in edge layers of tokamaks, where velocity shearing rate dv /dr approaches or even exceeds 
the drift wave frequencies, this theory has to be strongly modified. 

The experimentally relevant strong shear regime was studied analytically in numerous papers 
applying the gyrokinetic approach (see, for example, Refs. [10], [11]; we mention also Ref.[12] as an 
example of a recent work). In our papers (see Refs. [13]- [17]) the investigations of linear evolution 
of drift waves in shear flows applied Kelvin's method of shearing modes or a so-called non-modal 
approach. This treatment is performed in terms of a physical time variable without any prior spectral 
decomposition in time. Kelvin's method appears to be very effective for the analytical investigations 
of the temporal evolution of numerous low frequency modes in plasma shear flows[13]-[17]. Using 
one of the simplest models available describing three-dimensional plasma turbulence at the edge of 
a magnetic confinement device, the Hasegawa-Wakatani equations[18], we find that under condition 
of strong flow shear non-modal evolution of the drift- type perturbations supersedes the modal 
development of linear instabilities. Drift modes are suppressed before they can grow, preventing the 
development of nonlinear processes and imposing its own physical character on the dynamics. The 
suppression of the drift resistive instability in the case of sufficiently strong flow shear is a non-modal 
process, during which the initial amplitude of the separate spatial Fourier mode of the perturbed 
electrostatic potential decreases with time as (v' t)~ 2 . 

In this paper we extend the non-modal treatment of the drift modes stability onto the plasma 
flows with a spatially homogeneous time-dependent velocity shear. In Section II we consider linear 
and weak-nonlinear solutions of the nonlinear Hasegawa-Mima equation for plasma flow with strong 
time-dependent velocity shear. In Section III we discuss the correlation properties of drift waves 
in plasma shear flows and their influence on the saturation of the resistive drift instability. In that 
section we pay the attention to flows with a "modest" velocity shear (of the order of or greater 
than the instability growth rate, but less than the drift waves frequency). In this case the stage at 
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which nonlinear processes determine the saturation of the drift resistive instability occurs before the 
development of the linear non-modal suppression of drift instability. Conclusions are given in Section 
IV. 



II. Non-modal drift waves in plasma flow with strong time- 
dependent velocity shear 

We investigate the temporal evolution of drift modes in time-dependent shear flow using the 
Hasegawa-Wakatani equations [18]. We use Cartesian coordinates aligned so that z is in the di- 
rection of the mean magnetic field, y is in the direction of the mean flow, and x is in the direction of 
the inhomogeneity of the mean flow and plasma density, assuming that these coordinates correspond 
to toroidal (</?), poloidal (9), and the radial (r) coordinates, respectively, of the toroidal coordinate 
system. In this geometry the Hasegawa-Wakatani system of equations for the dimensionless density 
n = n/n e and potential = eip/T e perturbations (n e is the electron background density, T e is the 
electron temperature) is [18] 

9 ( 9 Tr , , d c ( dd> d d<b d \ \ ~ , d 2 . 

d T T . . d c ( d(b d d<h d\\ d<h d 2 . _ . 

in + Vo(x,t) — -- Uf— --^ — )n + v de ^ = a— (n - 0) , (2) 



dt dy B \dy dx dx dy J J dy dz 2 

Vq (x, t) is the velocity of the sheared flow, a = T e j n e 2 r] ll , r).. is the resistivity parallel to the homo- 
geneous magnetic field Byz, p s is the ion Larmor radius at electron temperature T e , Vd e = cT e /eBL n 
is the diamagnetic drift velocity, L^ 1 = — rflnn 0e (x) /dx. We transform Eqs.(l), (2) to new spatial 
variables £,77, 

t — t, £ = x, 7] — y — J V (x,t) dt, z = z. (3) 

which are the generalizations of the convective coordinates, used previously in systems with homoge- 
neously sheared stationary flows [13]. In these coordinates the linear convection terms in the system 
of equations (1), (2) are excluded, but the Laplacian operator A now becomes time-dependent, 

» £ id f dv (i,t) d\(d / avjtf.t) o\ » (4) 

dx* + dy'-> \d( J d( dri)\dZ J d( d, 1 ) + drf { > 

leaving us with an initial value problem to solve. This time dependence is responsible for the shearing 
of the waves pattern by the basic flow in the convective frame of reference. 

The analysis of the effect of the flow shear spatial non-uniformity on the temporal evolution 
of the electrostatic potential of drift modes has shown [17], that for the velocity profiles without 
inflection point this effect is subdominant at all stages of the drift modes evolution. Therefore we 
consider here the case of the spatially uniform time-dependent velocity shear, for which 

— — dt = v a (t) , (5) 
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where v' Q = const is a parameter with dimension of the velocity shear, a (t) is a function of time. 
Introducing Fourier expansion with respect to £, 77 and ^-coordinates with k±, I and k z as the 
wave numbers conjugate tothe spatial coordinates £, i] and z, respectively, the linearized system (1), 
(2) with omitted nonlinear convective terms may be combined into the following equation for the 
potential (t, k±, /, k z ): 

1 1 d 2 o 2 d 1 86 

In Eq.(4) C = ak 2 /p 2 l 2 v' = T e k 2 /p 2 l 2 v' n e 2 rj ll ^> 1, S — Iv^/vq. The Laplacian operator A is equal 
to 

A(r,k ± ,l) = -(k ± -lv' a(t)) 2 -l 2 . (7) 

In variables £ and rj the linear solution to Eq.(6) has a form 

(£, »7,*) = y dk ± J dl<t> (k ± , I, 0) g (k ± , I, t) e ik ^ +il \ (8) 

where 0(/cj_,/,O) is the initial data and g(k±,l,t) is the linearly unstable solution. It was shown 
in Ref. [13] (see, also [17]) for the case of the spatially homogeneous time-independent velocity shear, 
that for v' < 7, where 7 is growth rate of the resistive drift instability in plasma without shear flow, 
the solution g (k±, /, t) in times t < {y' Q lp s )~ l has an ordinary modal form, 

g{k^l,t)=e^ t ^ t (9) 

with known frequency of drift wave Retu = ood = Sv' / (1 + p 2 {k\ + l 2 )) , and growth rate of the 
resistive drift instability 7 = Imw = uj d p 2 {k\ + I 2 ) /ak 2 (1 + p 2 [k\ + I 2 )) ~ v' S 2 /C, where C = 
Cj {k\ + l 2 ). Using solution (9), which has a form as in plasma without any shear flow, weak 
nonlinear theory, which is grounded on the Hasegawa-Wakatani system (1), (2) or Hasegawa-Mima 
equation[19] may be easily developed [20]. This theory will be free from the known problem of the 
solutions singularities at critical level points x = x , where phase velocity of wave becomes equal to 
flow velocity, i.e. u — k y vo (xo) = 0. They emerge, when linear, as well as weak-nonlinear, solutions are 
obtained by using spectral transform in time and variables x, y of the laboratory frame of reference. 
It is worth to note, that it is unreasonable from the very beginning to apply the spectral expansion 
in time in the case of the time dependent flow velocity. 

In the laboratory frame solution (8) becomes nonseparable in space and time and therefore 
quite different from the normal mode assumption. In the case of the time dependent velocity shear 
it has a form 

(r, t) = J dk ± J dl<j) (k ± , I, 0) g (k ± , I, t) e i^-lv' a(t))x+ily _ ( 1Q ) 

Only in times, for which \v' a (t)\ 1 (in times t \vq\~ 1 in the case of stationary velocity shear 
) solution (10) has a normal mode form. Because of the time dependence k x = kj_ — lv' a (t) of the 
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wave number component along the flow shear, the modes in the laboratory framebecome increasingly 
one-dimensional zonal-like as the perturbed E x B velocity tilts more and more closely parallel to 
y-axis. 

The modal solution (8), (9) follows over times a(t)v' lp s > 1 strong non-modal evolution of 
<f>(t,k±,l), resulted from the enhanced by shear flow time-dependent dispersion, which does not 
allow the weak nonlinear three waves decay processes. For strong flow shear case, when lvd e % v' , 
the first term of Eq.(6) over times t > (lvd e )~ l may be omitted[17]. In that case Eq.(6) reduces to 
the Hasegawa-Mima equation[19], which with convective nonlinearity (omitted in Eq.(6)) in the case 
of the time- dependent velocity shear has a form 

[a (k±,l,t)(f)(t, k±,l)\ +ilv de <f)(T, k ± ,l)=p 2 s ^ J dk ±1 J dh J dk ±2 J dl 2 A (t, k 2± ,l 2 ) 

x (k 1± l 2 - k 2± h) (j) (t, k 1± , h) (t, k 2± , l 2 ) 6 (k± - k 1± - k 2± ) 6(1 — h — l 2 ) , (11) 



where A (r, k 2 ±, l 2 ) is determined by Eq.(7)) and 

a (k ± , l,t) = l + p 2 s (l 2 + (k ± - v' a (t) if^j . 



(12) 



The approximate weak nonlinear solution to Eq.(ll)) for the potential is easily obtained in the 
form of the power series in the amplitude of the initial data </>(£, Z,0), <f>(r,k±,l) = <po(T,k±,l) + 
0i ( T ,k_L,l) + where O and 0i are linear and square terms with respect to (£, I, 0), with 



0o (k ± , /, t) = (k ± , I, 0) g (k ± , I, t) , 



(13) 



and 



0i (k±,l,t) 



Ba(k±,l,t) 



dtiexp 



to 



tl 

—ilvde J dt 2 a~ x (k±,l,t 2 ) J dk±i J dl\ J dk± 2 J dl 2 



to 



x A (ti, k 2± , l 2 ) (k 1± l 2 - k 2 ±h) (k 1± , h, 0) (k 2± , l 2 , 0) g (k 1± , h,t) g (k 2± , l 2) t) 

x6(k ± - k 1± - k 2± ) 6(1 -h- l 2 ) , 



(14) 



where 



g(k ± ,l,t) 



a(k±,l,t) 



cxp 



—ilv, 



de 



t 

J dUa' 1 



(k±,l,h) 



to 



(15) 



It follows from Eqs.(13) and (14) that over times, for which \a(t)v' l\ p s > 1 the initial per- 
turbations of the potential with given lp s will be suppressed before the development of the resistive 
drift instability, when non- modal evolution begins at times t < 7 _1 . For the particular case of the 
time-independent flow shear solution (10) was obtained and discussed in Ref.[13]. For power-like 
time dependence of the velocity shear with <9Vo/<9£ = v' a n t n , a (t) = a n t n+1 /(n + 1) and 

*(k ± .l.t)= i p: I I' 2 ■ (k. -r,V^i^-VY (10) 



(n + 1) 



the non-modal development begins at times t > (v' a n lp s ) at which linear O (t, k±, I) and weakly 
nonlinear <pi (t, k±, I) parts of the potential rapidly decay with time. The fluctuation spectrum for 
long time will be dominated by Fourier components with sufficiently small values of wave number 
component / for which condition \a(t)v' l \ p s > 1 is not met during time considered. 

III.Renormalized hydrodynamic theory for drift modes in 
plasma shear flow 

With variables £,77, Eqs. (1) and (2) have a form 

*(i-i(g£- IS) **-•£<»-«• < i7 » 

d c ( d(j) d d(j) d \\ , d(f) d 2 

dt-B\W v Wr^)) n + Vde ^y =a d^ {n ~^- ( 8) 

The above calculations show that with variables £ and 77 we exclude from Eqs.(l) and (2) the spatial 
inhomogeneity originated from shear flow. That gives an opportunity to investigate in linear approx- 
imation the temporal evolution of the separate spatial Fourier mode of the electrostatic potential 
with definite wave numbers k± and /. It is interesting to note that transformation (3) conserves 
the E x B convective nonlinear derivative in Eqs.(l) and (2) in the form similar to one in a plasma 
without any flows. The flow shear contains only in time dependent Laplacian operator and become 
apparent in linear non-modal effect of the enhanced dispersion during times \a(t)v' l\ p s > 1 ( or, 
at the case of stationary velocity shear, when t > (fg/ps)" 1 ) after the stage of the ordinary modal 
evolution[13]. With new variables £1,771, which are determined by the nonlinear relations 

to to 

the convective nonlinearity in Eqs. (17) and (18) becomes of the higher order with respect to the 
potential 0, 

c / d<fi d 90 9 

B \ dr] <9£ <9£ <97/ 

90 (t) f d 2 0(*i) d<p(t) 




dt! 



<9t7<9£ 

to 



<97/<9£ 
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<9£ 2 J drf \ J drjdS, 



Omitting such nonlinearity, as well as small nonlinearity of the second order in the Laplacian, resulted 
from the transformation to nonlinearly determined variables £1,771, we come to linear equation (6) 
with solution (8), where wave numbers k±,l are conjugate there to coordinates £1, 77! respectively. 
With variables £ and 77 this solution has a form 

0(£,77,O = Jdk ± J di0(A; ± ,Z,O) 5 (A; ± ,Z,ii)e^ + ^-^>)-*), (21) 

where 

to *o 

Eq.(21) is in fact a nonlinear integral equation for potential 0, in which the effect of the total 
Fourier spectrum on any separate Fourier harmonic is accounted for. The functions £ (t) and 77 (t) 
in the exponential of Eq.(21) involve through Eq.(22) integrals of 0, which in turn, involve in their 
exponentials the integrals (22) and so on. This form of solution, however, appears very useful for the 
analysis of the correlation properties of the nonlinear solutions to Hasegawa-Wakatani system and 
for the development of the approximate renormalized solutions to Hasegawa-Wakatani system, which 
accounted for the effect of the turbulent motions of plasma on the saturation of the drift-resistive 
instability. Consider first the dispersion tensor of random displacements £ (t) and 77 (t) of the plasma 
resulted from its motion in such turbulence. The variances of these displacements are determined by 
the relations 

k«m = k«w = (((&)') = £f*>)*{^^)> P3) 

to to 

Kin (t , k) = K(n {t) = ( (mm) ) = 4 /*,/*, ( gffi °*M ) , (24 ) 

to to 



K m (Mo) = K m (t) = ((77 (t)) 2 ) = ^Jdt 1 J dt 2 



c 2 L L /90(*iW(t2) 



<9£ <9£ 



(25) 



to to 

where brackets (...) mean an average over ensemble of random values of (k±, I, 0). Using in Eqs.(23)- 
(25) solution (21) and the Corrsin's[21] independence hypothesis, 



x 



(h±, Ml) (h ± , l 2 , t 2 ) jk^+Hri-ik^-ilW^ = ( ( fcl _ L> kj 0) (fc 2± , Z 2 , 0)) 

9 (*i±, Zi, tx) (fc 2± , Z 2 , t 2 ) / e «i-L(«*i)-«fa))+«i S (k 1± + fc 2± ) 5 (h + Z 2 ) , (26) 



we find for the components (t) , K^ v (t) and K m (t) of the dispersion tensor the expressions 



2 

K tt ^ = hj dtl / dh J dk± / d l\(l>(k±,l,0)\ 2 l 2 9(k±,l,t 1 )g(-k ± ,-l,t 2 ) 



to to 



x ( e 



= ifc 1± (?(tl)-5(t 2 ))+i«l(??(tl)-5?(t2)) 



(27) 



t t 

K^{t) = J dtx J dt 2 J dk ± J dl\(f)(k ± ,l,0)\ 2 k±lg(k±,l,ti) g (-k ± ,-l,t 2 



t() to 



x ( e 



3 ifcix(£(ti)-£(*2))+iMi7(ti)-i7(t2)) 



(28) 



2 

K vv( t ) = ^ J d h J dt 2 J dk ± J dl\<P(k ± ,l,0)\ 2 l 2 g(k ± ,l,t 1 )g(-k ± ,-l,t 2 ) 



to to 



x < e 



3 ifcix(£(ti)-£(*2))+iMi7(ti)-i7(t2)) 



(29) 



During times |a(t)fo/| p s < 1, the solution g to the Wasegawa-Wakatani system is of a modal form 
(9). Assuming that the displacements £(t), obey the Gaussian statistics with mean zero, we 
find in this case the following relation for (t): 



t t 



K tt ^ = ^2 J dtl J dh / dk± / dl\<j){k ± ,l,0)\ 2 l 2 exp ( 7 (t 1 +t 2 )+iou d (h - t 2 )) 



to to 

x exp 



-k 2 ± Kx (h, t 2 ) - kJKtr, (h,t 2 ) - h 2 K m (h,t 2 ) 



(30) 



where 



Ktt(t 1 ,t 2 ) = (tt(t 1 )-t(t 2 ) 



Ks v (t u t 2 ) = ( (£ (h) - £ (t 2 )) ft fa) - v(t 2 )) } , 
K m (t 1 ,t 2 ) = ((?](t 1 )-?j(t 2 )) 2 ). 



(31) 
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With time variables r = t\ — t 2 and t — (tx + t 2 ) /2, Eq.(30) becomes 
o t+| t t~i 



(t) = 



B 2 



dr dt+ dr dt 



dk ± J dl\(j)(k ± ,l,0)\ 2 l 2 exp (2~f(k±,l)t + iu d (k±,l)T) 



x exp 



2B 2 



J dn J dh J dk 1± j dli |0 h, 0)| 2 |[k_L x k 1± ]| 2 



t-|(r+n) 



x exp ( 27 (feu., Zi ) ti + (feu., h ) Ti - ^ (k 1± • r ^ + ^ - k 1± • r (t^ - ^ 

?+i(r-n) 

y dn J dh J dk 1± J rf/i|0(A; 1± ,/ 1 ,O)| 2 |[k ± xk 1± ]| 2 



2B 2 



f-i(r-n) 



x oxp I 27 Zi ) ti + (fcu.,Zi ) n - ^ ( (ki_L • r (^i + -| 



- ki ± • r (ix - 



Tl 



(32) 



where r = (£,77) and kj_ = (k±,l). The similar relations (see Eqs.(27)-(29) above) may be obtained 
for K m and K^ v , changing I 2 with k\ for K m and with — k±l in the integrand for K^ v , respectively. 
Eq.(32) contains two time scales. The fast time scale, determined by lu^ 1 , is the correlation time 
and the slow time scale, determined by 7 -1 may be considered as the relaxation time. In the case, 
when 7- 1 > a;" 1 , Eq.(32) may be substantially simplified. As it follows from (32), the integrals over 
slow time t\ and infinite sequence of integrals, which appear in r (t\ ± y) over slow time scales, are 
calculated over narrow time intervals of the order of r ~ T corre i at i on ~ uo^ 1 over which the variation 
of the integrands that occurs from the slow time is negligible. In such case the integrals over t\ may 
be approximately calculated as 

t+i(r±n) 

f d^/(Qw(r±n)/(*). (33) 



t-§(T±n) 

With quasi-Markovian approximation, 



\ ((k 1± -r (h + ^) - k 1± • r (£ - |)) 2 \ « nC7 



(34) 



which, in fact, resulted from Eq.(33), we obtain the integral equation for the function C (k±,l,T): 

X - {k\K^ (*!, t 2 ) + kJKfr (*i,t 2 ) + / 2 ^ w (ti,t 2 )) 

= 5 (*i^€ (* + + fe±^€, (* + I) + l 2 K m (t+ T -,t- I)) 

= C (Jfej., Z,t) r 

00 

= r-J I dh\<p(k 1± ,h,t)\ 2 \[k^xk 1± ]\ 2 J ^eN^iJiKM^n 



C 2 /" f 1 , -M2 lr ll2 C Zi,^) 

~ T ~^ dk i± dhUlk^l^t)] kiXkn — ^- — f, 



(35) 



where |0 h, t 



(k 1± , h, )| 2 e 2 ^ kl± ' h $. Then K# (t) becomes equal to 



c 

~B 2 



o t+| t *-§ 
/ dt + I dr I dt 



x exp 



y dfcn. y dZi|0(fci ± ,Zi,O)| 2 e^ fcl - L ' ,1 ) f |[k J .xk 1± ]| : 



dA; ± / rf/|0(A; ± ,/,O)| 2 / 2 e 27{ ^' i) * + ^ ( ^' /)T 



C(k 1± ,h,t) 



2c 2 /" [ 2 f ~\ / 

■jp / c/A:j_ / dll / (it |0 /, t 



2 C(k±,l,t) 
uj 2 d (k ± , I) 



(36) 



With Eq.(36) and similar equations for K^ v (t) and K m (t) we obtain a more general equation, 

klK^ (t) + 2kJK iv (t) + l 2 K m (t) 
t ^ t 

= ^J dk 1± J dh J dt\<t> (h ± ,h,t)\ 2 \[k ± x k 1± ]| 2 ^^fj = 2 J dtC(k ± ,l,t)- (37) 
o d ' o 

With Eq.(37) we obtain the renormalized form of the potential (21), in which the average effect of 
the random convection is accounted for, 



<t> (£,r],t) = / dk± / dl(j) (k±, 1, 0) e 



odt+yt-f dtC(k±,l,t)+ik±Z+ilr) 
o 

The saturation of the instability occurs when <9 (0 (£, 77, t)) 2 /dt = 0, i.e. when 

7(k±,l) = C(k ± ,l,t) 

,2C(k 1± ,h,t) 



(38) 



y dfcxj. y ^1 10 /i,t)i 2 i[k ± x ku_]i : 



(39) 



u\ (ki±,h) 

From the double Eq.(39) we obtain the equation, which determines the level of the instability satu- 



ration 



7(*l,0 „ 2 



c 



y dk 1± J dh\(f) (k 1± , h, t) | 2 1 [k ± x k 1± ] | 



2 7(^i±^i) 



(40) 



Wd(fci±,Zi)" 

The sought-for value in Eq.(40) is a time t sat at which the balance of the linear growth and nonlinear 
damping occurs for given initial disturbance (ki±, h, 0) and dispersion. With obtained t sat the 
saturation level will be equal to \<p{t sat )\ 2 ~ / dfcu. / dZi |0 (k 1± , l u 0)| 2 e 27 ^ Jl )*«*«. Also, the well 
known order of value estimate [20] for the potential in the saturation state is obtained easily from 
Eq.(40), 



e0 



T e k±L n 



(41) 



Obtained results show that the nonlinearity of the Hasegawa-Wakatani system of equations 
in variables £ and 77, with which frequency and growth rate are determined without spatially inho- 
mogeneous Doppler shift and wave number is time independent, does not display any effects of the 
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enhanced decorrelations provided by flow shear. In the laboratory frame of reference such spatial 
Fourier modes are observed as a sheared modes with time dependent component of the wave number 
k x = k± — lv' a (t) directed along the velocity shear, 

(r, t) = J dk ± J dl<f> (k ± , I, 0) e *(*x-K«W)*-Ww+^d*+7t-ifcx«ti)-^(ti)_ (42) 

To make the next analysis simpler we consider in what follows the case of the spatially homogeneous 
time-independent velocity shear, a (t) = t. The displacements £ (t) and rj (t) are observed in the 
laboratory frame as the displacements x (t) and y (t) which are equal to 

t 

x(t)=m = ~J^dt 1 



dy 

to 



t 



= -i^ J dh J dk ± J dl(P(k ± J,0)le^ dtl+ ^ tl+i{k± - l < tl)x+ily ~^ (43) 



to 



and 



t t 



y(t) = Jv~(t 1 )dt 1 = ^Jdt 1 ^ + jdh-^^xih) 

to to *0 

t 

= J dk ± J dl J dt!(f)(k ± ,l,0) (k±-lv' (f-fi)) 

to 

xe jaJ d tl+7tl+ i ( fc ±- to (')tl) :E +^?/-* fe x5(tl)-^»?(tl) (44) 



Applying the procedure presented above to the calculations the variants of the plasma dis- 
placements K xx (ti,t 2 ) in variables x, y we find that K xx (ti,t 2 ) has the same form as (ti,t 2 ) 
in Eq.(36) in which, however, Ud is changed by Ud — lv' x, assuming that ujd — lv' x ~ Ud, and 
C (k\±, li, t) is changed by C\ (k±±, l\, t,x), determined by Eq.(39), in which the replacement ujd by 
^d — lv' Q x ~ Ud has to be done. The correlation K yy (t) is 



c 2 



t+§ t t-§ 



J dr J dt + J dr J dt J dk ± J dl \<j> (k ± , I, t) \ 
x (k ± - lv' (t -t- T ^j) (k ± - lv' (t - t + 



2 i(w d (k±,l)-lv' x)T-Ci\T\ 



* 2 ^Jdk ± JdlJ dt\d> (k ± , l,t) I 2 (k ± - lv' (t - t)f ^ )2 . (45) 

The calculation of the integral over t needs the knowledge of the C\ (t) dependence. In the case of 
the stationary perturbations with 7 = the integral over t is easily calculated and we obtain for 
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times t > T corre i ation ~ (uj d 



lv' Q x) 1 the well known result [22] 



K yy (t) 




Ci (k±,l,t,x) 



(u d (kj_,l) - lv' x) 2 



(46) 



which displays the effect of the anisotropic dispersion conditioned by flow shear, observed in the 
laboratory frame of reference - dispersion increases much faster along flow than in the direction of 
the flow shear. 



In this paper we present the results of the analytical investigations of the linear and nonlinear 
evolution of drift modes in shear flows with time dependent velocity shear. The transformation (3) 
to convective frame of reference gives an opportunity to obtain the exact solution for time dependent 
flow shear. The investigations demonstrate that there are different time scales in the dynamics of drift 
waves in shear flows. Right with convective variables £, r\ (3) during times, for which \a (t) v' l \ p s <C 1 
(or, for the stationary velocity shear, when t <C (Ipsv'o)' 1 ) the solution (9) of the modal type is 
obtained. That modal solution already during early times, for which |a:(£)i>o| > 1, is observed in 
the laboratory frame as a sheared mode with time-dependent wave number and frequency. In this 
time domain in the convective variables the dispersive properties of the plasma displacements in 
drift turbulence are the same as in plasmas without shear flow. Using the developed two-time scale 
procedure of the calculation of the dispersion tensor of random displacements of the plasma we obtain 
the renormalized form (38) of the solution for the electrostatic potential of the unstable drift waves. 
This solution gives the nonlinear integral balance equation (39), which determines the level of drift 
turbulence in the steady state, resulted from random turbulent motion of plasma in the unstable drift 
turbulence. Balance equation (39) governs the process and level of the drift turbulence saturation. 
It does not include any effects associated with flow shear over times during which the solution (9) is 
valid. 

The effect of the anisotropic dispersion of plasma displacements is presented by Eq.(46). It 
presents the variance of the plasma displacements resulted from the modes of the ordinary modal 
form (9), observed in the laboratory frame. This effect has nothing in common with effect of the 
"enhanced suppression of the instability by shear flow". The saturation of the instability, as it was 
demonstrated above, resulted from the balance (39) of the growth rate and nonlinear damping in the 
frame of reference, where the solution in the ordinary modal form is determined. 

This result may be extended to other fluid models of plasma. The fluid equations, obtained 
in drift approximation, in which all nonlinearities other than E x B nonlinearity are ignored, do not 
include among the nonlinear mechanisms of the instability saturation the process of the enhanced 



VI. Conclusions 
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nonlinear decorrelation by velocity shear. It is important to note, that the same conclusions con- 
cerning "suppression of the instability by the enhanced nonlinear decorrelation by velocity shear" are 
completely applicable to stratified shear flows of incompressible fluids, where the nonlinear convective 
derivative 

dip d dip d , . 

^ dz dx dx dz 

presented in terms of the stream function ip = e y ip, with which v = y x ip and where basic flow 
along axis x with velocity gradient along axis z is assumed. Eq.(47) also conserves its form under 
transformation (3). 

The process of the instabilities suppression by the enhanced nonlinear decorrelation by velocity 
shear is treated appropriately with nonlinear kinetic (gyrokinetic) theory of plasma with shear flow 
(see, for example, nonlinear gyrokinetic turbulence simulations of Ex B shear quenching of transport 
in Refs.[6], [7]; we mention also Ref. [23], [24, 25] as examples of a recent works). We have obtained 
analytically in Ref. [25] that nonlinear effect of the suppression instability by shear flow deduced 
straightforwardly in Vlasov-Poisson system of equations for the magnetized plasma as the effect of 
the enhanced by shear flow nonlinear resonance broadening and it may be dominant in the processes 
of the instabilities saturation. 
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